Relaxation oscillations in long-pulsed random lasers 
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We have measured the evolution of the light intensity of a random laser during a nanosecond 
pump pulse. Relaxation oscillations in a titania random laser were observed in the time trace of the 
total emitted intensity. We compare our experimental results with a simple model, based on the 
four-level rate equations for a single mode laser. 

PACS numbers: 42.60.Rn, 42.25.Dd, 42.55.Zz 
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I. INTRODUCTION 

Relaxation oscillations of conventional lasers are a well- 
understood phenomenon. [l| They are especially impor- 
tant for continuous wave and long-pulsed lasers. In a ran- 
dom laser, a medium in which gain is combined with mul- 
tiple scattering of light, relaxation oscillations can also 
occur, as was predicted by Letokhov in 1968. [2j A pre- 
eminent experimental demonstration of a random laser 
was published by Lawandy et al. in 1994 followed 
by many others [4|, [5j, Is III ■ Only recently Soukoulis and 
coauthors have presented measurements of relaxation os- 
cillations in single modes of a picosecond pumped random 
laser system. [8[ Different numerical calculations on ran- 
dom lasers also show this oscillatory behavior. [9T [ToL fliT] 
To our knowledge, no measurements have been per- 
formed on relaxation oscillations in random lasers in the 
interesting regime of long pulses. 

In this report we present measurements of the time 
evolution of a nanosecond pumped random laser system. 
We compare our experimental observations with a simple 
model, based on the four-level rate equations for a single- 
mode laser. 



II. A SIMPLE MODEL 

The evolution of the excited molecules Ni and the 
number of photons q in a laser are described by the well- 
known four-level rate equations [l[ 
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where P is the pumping fluence that is absorbed by the 
molecules inside the cavity, t c the cavity decay time, r 
the spontaneous-emission life time, and (3 the beta factor 



of the laser. The beta factor is defined as the amount 
of spontaneous emission that contributes to the lasing 
mode, and can also be determined for a random laser [l2T ] . 

From these rate equations Woerdman and coauthors 
have derived an equation for the frequency of the re- 
laxation oscillations explicitly including the spontaneous 
emission: 1 13] 
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where o; res is the relaxation oscillations frequency, and 
M the scaled pump fluence, defined as the ratio of the 
absorbed fluence P and the threshold fluence Pth- 

We apply this model to the multi-mode pulsed random 
laser by simply changing (3 to /3 mm and r c to r c mm , i.e. 
we use a mean cavity decay time and a mean beta factor 
to describe our multi-mode random laser. This simpli- 
fication accurately describes the threshold behavior of 
random lasers [14| . 



III. EXPERIMENTAL APPARATUS 

The random laser used in our experiments consists of a 
suspension of TiC>2 particles (mean diameter of 180 nm, 
volume fraction of 10%) in a solution of Sulforhodamine 
B in methanol (1 mmol/liter; pump absorption length, 
104 /im; minimal gain length, 83 /im [HI). The suspen- 
sion is contained in a fused silica capillary tube, with 
internal dimensions 100 x 2 x 2 mm 3 . To measure the 
mean free path of light in this sample, we performed an 
enhanced-backscatter cone experiment [l6j and an escape 
function experiment [l7l |. We found a transport mean 
free path of 0.46 ±0.1 ^m at 633 nm (effective refractive 
index, 1.48±0.04). 

The samples were excited by a pump pulse at 532 nm, 
provided by an optical parametric oscillator (OPO) 
pumped by a Q-switched Nd:YAG laser (Coherent In- 
finity 40-100/XPO). The pump pulse had a duration 
of 3 ns and a repetition rate of 50 Hz. The pump 
light was focused with a microscope objective (water- 
immersed, numerical aperture NA, 1.2) onto the sample 
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(focus area, 12 ±6 /im 2 ), reaching an intensity in the or- 
der of 1 mJ/mm 2 . The central wavelength of the emitted 
light is 594 nm [3| , and the narrowing factor (defined as 
the spectral width of the emitted light far above thresh- 
old divided by the spectral width far below threshold) 
is 8. The light emitted by the random laser was col- 
lected by the same microscope objective. The emitted 
light was detected by a 25 GHz photodiode (New Focus 
1404), read out by an oscilloscope (Tektronix TDS 7404, 
analog bandwidth, 4 GHz). The resulting time resolution 
was 100 ps. To obtain a good signal-to-noise ratio, the 
data shown is an average of 100 oversampled time traces. 
The pump light was filtered out of the detection path 
by use of a colored glass filter with an optical density of 
more than 4 at the wavelength of the pump laser. 



IV. MEASURED RELAXATION 
OSCILLATIONS 

The normalized time trace of the pump pulse and the 
normalized time trace of the total emitted light from the 
random laser far above threshold are shown in Fig. [T] 
Overall, the duration of the pump laser pulse is longer 
than the duration of the pulse of light the random laser 
emits. We see in the pulse emitted by the random laser 
first a fast decay, followed by a slower exponential de- 
cay. The fast decay is due to the stimulated emission 
in the random laser. In the second part of the decay 
the population inversion is no longer present, and the 
spontaneous emission causes a slower decay of intensity. 
These observations are in agreement with other random 
laser experiments. @ In Fig. [1] relaxation oscillations in 
the emitted light are clearly visible near the peak inten- 
sity. 

We measured the time evolution for different input 
fluences. In Fig. [5] the normalized intensity is plotted 
versus time for four different pump fluences 18[. The 



time traces are shifted vertically with respect to each 
other for clarity. The time trace at a pump fluence of 
0.06 mJ/mm 2 is below threshold, while the time traces 
with higher pump fluences are above threshold. We ob- 
serve that relaxation oscillations occur above threshold 
and become more pronounced when the pump fluence 
increases. 

The frequency of the relaxation oscillations are com- 
puted from the time traces. We determine the times at 
which the intensity is at a local maximum. The differ- 
ence of two consecutive local maxima At is the period, 
and the frequency of the relaxation oscillation v rc \ is given 
byl/(Ai). 



V. COMPARISON OF THE MEASUREMENTS 
WITH THE MODEL 

We have inferred the relaxation oscillation frequency 
for different pump fluences from our measurements. In 
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FIG. 1: Measured time traces of the pump pulse (gray) 
and emission output above threshold (black, input fluence 
= 0.47 mJ/mm 2 ) of a titania random laser. The pump pulse 
duration is much longer than the duration of the emitted light 
of the random laser above threshold. Relaxation oscillations 
in the emitted light are clearly visible near the peak intensity. 
The decay time of the emitted light is first dominated by stim- 
ulated emission. In the second part of the decay-curve, the 
spontaneous emission is dominating. 



Fig. [3] the measured relaxation oscillations frequencies 
are plotted versus the scaled pump fluence M. The 
relaxation-oscillation frequency significantly decreases 
when the scaled pump fluence increases from 1 to 2. A 
further increase of the scaled pump fluence does not sig- 
nificantly change the frequency of the relaxation oscil- 
lations. The result of Eq. ^ is depicted for different 
cavity decay times. This cavity decay time is the only 
parameter that could not be directly determined by our 
experiment. The trend of the model for a fixed cavity 
decay time is that, in contrast to our measurements, the 
relaxation-oscillations frequency increases for increasing 
pump fluence. Only for large (> 3) normalized pump flu- 
ence the fit of the model for a cavity decay time of 5 ps 
fits reasonably. 

The difference in the observed trend of the relaxation- 
oscillation frequency between the simple model and our 
experiments, is probably due to the size of the gain vol- 
ume: The gain inside the gain volume of the random 
laser saturates at threshold [l4] . A further increase of 
the pump fluence will lead immediately to a larger gain 
volume, and a corresponding longer cavity decay time. 
Apparently, this behavior saturates above a scaled pump 
fluence of 3. 

We have used the distribution of the phase-delay time 
to determine the mean cavity decay time Tc, given by [3| 
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FIG. 2: Measured time traces of the emission intensities of 
a titania random laser for four different pump fluences. The 
traces are vertically shifted with respect to each other for 
clarity. Relaxation oscillations become more pronounced at 
higher pump fluences. 



The measured relaxation oscillations were compared 
with a simple model, based on a single-mode continuous- 
wave laser system. The observed trend of our experi- 
ments differs from the expected trend, due to the increase 
of the gain volume and the corresponding cavity decay 
time for an increase of the scaled pump fiucncc. The cav- 
ity decay time determined with the fit from the simple 




with L the length of the gain volume and, for 
non-resonant scattering, the diffusion constant D = 
Co£/(3n'), where cq is the speed of light in vacuum, £ 
the transport mean free path and n' is the real part of 
the effective refractive index of the medium. For our ti- 
tania random laser we find a 7c in the order of 10~ 13 s, in 
contrast to 10~ 12 s that is suggested by the agreement in 
Fig. [3] at M = 3. The difference between the two cavity 
decay times is a factor 10. This deviation could origi- 
nate in part from the difference between the mean cavity 
decay time of all modes and the mean cavity decay 
time of the lasing modes: Our experiment suggests that 
random lasing preferentially takes place in modes with a 
much longer than average cavity decay time. 

VI. CONCLUSIONS 



FIG. 3: Measured relaxation oscillations frequency as a func- 
tion of the normalized pump fluence (squares). The simple 
model [Eq. (J5J)] is plotted for different cavity decay times: 0.7 
(dash-dotted line), 2 ps (dotted line), and 5 ps (dashed line). 
The fit of the model for a cavity decay time of 5 ps fits reason- 
ably for high pump fluences, but corresponds to a surprisingly 
large value of r. 

model is a factor 10 higher than the mean cavity decay 
time of our sample. Our experiment suggests that the 
modes contributing to random laser emission have a cav- 
ity decay time much longer than the average cavity decay 
time. 



We have seen relaxation oscillations in our random 
laser, while looking at the time evolution of the total 
emitted light for different realizations of the sample. Mul- 
tiple modes contributed to these time traces, and we av- 
eraged the time traces over several realizations of disorder 
of our random laser sample. The resulting time trace still 
showed relaxation oscillations: a weighted average of the 
oscillations of all the underlying modes. 



Acknowledgements 

This work is part of the research program of the Sticht- 
ing voor Fundamenteel Onderzoek der Materie (FOM), 
which is financially supported by the Nederlandse Or- 
ganisatie voor Wetenschappelijk Onderzoek (NWO). We 
thank Boris Bret for the work that lead to the estimation 
of the transport mean free path. 



[1] A. E. Siegman, Lasers (Stanford University, University [3] N. M. Lawandy, R. M. Balachandran, A. S. L. Gomes, 

Science Books, Mill Valley, CA, Oxford, 1986). and E. Sauvain, Nature (London) 368, 436 (1994). 

[2] V. S. Letokhov, Sov. Phys. JETP 26, 835 (1968). [4] D. Sharma, H. Ramachandran, and N. Kumar, Opt. Lett. 



4 



31, 1806 (2006). 
[5] A. Rose, Z. Zhu, C. F. Madigan, T. M. Swager, and V. 

Bulovic, Nature (London) 434, 876 (2005). 
[6] H. Cao, Y. G. Zhao, H. C. Ong, S. T. Ho, J. Y. Dai, 

J. Y. Wu, and R. P. H. Chang, Appl. Phys. Lett. 73, 

3656 (1998). 

[7] S. Mujumdar, M. Ricci, R. Torre, and D. S. Wiersma, 

Phys. Rev. Lett. 93, 053903 (2004). 
[8] C. M. Soukoulis, X. Jiang, J. Y. Xu, and H. Cao, Phys. 

Rev. B 65, 041103 (2002). 
[9] D. S. Wiersma and A. Lagendijk, Phys. Rev. E 54, 4256 

(1996). 

[10] G. A. Berger, M. Kempe, and A. Z. Genack, Phys. Rev. 

E 56, 6118 (1997). 
[11] M. A. Noginov, J. Novak, and S. Williams, Phys. Rev. A 

70, 063810 (2004). 
[12] G. v. Soest and A. Lagendijk, Phys. Rev. E 65, 047601 



(2002). 

[13] J. P. Woerdman, M. P. van Exter, and N. J. van Druten, 

Adv. At., MoL, Opt. Phys. 47, 205 (2001). 
[14] K. L. van der Molen, A. P. Mosk, and A. Lagendijk, Phys. 

Rev. A 74, 053808 (2006). 
[15] G. van Soest, Ph.D. thesis, University of Amsterdam, 

www.randomlasers.com, 2001. 
[16] M. P. van Albada and A. Lagendijk, Phys. Rev. Lett. 55, 

2692 (1985). 

[17] J. Gomez Rivas, D. H. Dau, A. Imhof, R. Sprik, B. P. J. 
Bret, P. M. Johnson, T. W. Hijmans, and A. Lagendijk, 
Opt. Commun. 220, 17 (2003). 

[18] A slight time shift is possible, since the triggering of the 
oscilloscope was found to depend on the pump laser in- 
tensity. 



